Noise-Driven Mechanism for Pattern Formation 
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We extend the mechanism for noise-induced phase transitions proposed by Ibahes et al. [Phys. 
Rev. Lett. 87, 020601-1 (2001)] to pattern formation phenomena. In contrast with known mech- 
anisms for pure noise-induced pattern formation, this mechanism is not driven by a short-time 
instability amplified by collective effects. The phenomenon is analyzed by means of a modulated 
mean field approximation and numerical simulations. 
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I. INTRODUCTION 



Starting with the seminal work of Horsthemke and 
Mansour on the Verhulst model [|l], noise-induced phe- 
nomena have been a subject of intense interest g. Much 
of the early work dealt with noise-induced phenomena 
in zero-dimensional systems. More recently, it has be- 
come widely recognized that the effects of fluctuations 
on systems with a large number of degrees of freedom, 
so-called spatially extended systems consisting of coupled 
zero-dimensional systems, can be even more striking [[|. 
First-order ||, || and second-order || [?], || noise-induced 
phase transitions, noise-induced patterns || |l(| O, p"2| , 
and doubly stochastic resonance |D| are examples that 
illustrate the broad interest in the subject. The noise 
intensity is the parameter that controls the creation of 
order in all these processes. The underlying mechanism 
that drives the order in these extended systems is a noise- 
induced short-time instability that is amplified by collec- 
tive effects that do not occur in the absence of coupling 
or in the absence of noise. Moreover, as the intensity 
of the fluctuations increases, a re-entrance phenomenon 
occurs: the noise that first drove the system to an or- 
dered state now restores the disorder. It is important to 
recall that this particular mechanism is in a sense op- 
posite to that responsible for noise-induced phenomena 
in zero-dimensional systems [|[ |7| in that these noise- 
induced phase transitions are observed only when the 
associated zero-dimensional units exhibit no interesting 
noise-induced behavior. It was therefore thought for a 
long time that coupling of zero-dimensional units that 
undergo a noise-induced transition would exhibit no in- 
teresting collective effects. 

However, Ibahes et al. (T^| introduced a class of exactly 
solvable models that exhibit both noise-induced tran- 
sitions in the zero-dimensional case and noise-induced 
phase transitions in the associated extended system. 
They stress that this phase transition arises from an ef- 
fective equilibrium potential in the steady state and does 
not require a short-time instability or any other reference 
to the short time behavior of the system. The key ingre- 
dient is the combination of relaxations that rely on field 
dependent kinetic coefficients and the disordering effects 



of external fluctuations. Furthermore, in this case no re- 
entrance phenomenon occurs: the system becomes more 
ordered with increasing noise intensity. 

In this paper we extend this mechanism to pattern 
formation. We modify the coupling term of the model 
in such a way that a morphological instability can ap- 
pear. We will show, by means of numerical simulations 
and analytic calculations, that an increasingly ordered 
spatial structure develops as the intensity of the fluc- 
tuation increases beyond a critical value. Moreover, we 
show that for sufficiently high noise intensity there is non- 
monotonic behavior as a function of the coupling, that is, 
there exists and optimal value of the coupling for which 
the ordered structure exhibits maximum coherence. 

The paper is organized as follows. We introduce the 
model in Sec. [n] and explicitly show that no short-time 
instability can drive pattern formation in this system. In 
Sec. Ill and Appendix ^ we present a modulated mean 
field theory and the phase diagram of the model obtained 
from this theory. Numerical simulations that confirm the 
qualitative validity of the theoretical results are presented 
in Sec. Jtyj. Finally, we summarize our main conclusions 
in Sec. rV| 



II. THE MODEL 

As introduced in (lij, we consider the following 
Langevin dynamics for a space and time dependent scalar 
field <j) r (i): 

4>r =T{4> r )[-a4> r +C4> r ] + [T{dp r )] 1/2 i(r,t). (1) 

Here the local force — a<\> arises from a monostable local 
potential a(j) 2 /2 and £ stands for the spatial coupling 
operator. The field-dependent kinetic coefficient is 



i 



1 + C(j) 2 



(2) 



The coupling of this kinetic coefficient with the noise fa- 
vors fluctuations in the disordered state = 0. Both 
a and c are positive constants. The noise term is as- 
sumed to be Gaussian with zero mean value and correla- 



2 



tion function 



<£(r,t)£(r',0> =2a 2 S(t-t')S(r 



(3) 



where the brackets ( ) denote a statistical average. 
The associated stochastic integral is interpreted in the 
Stratonovich sense 0. 

In the absence of coupling (£ = 0), the model (|l|) re- 
duces to a collection of zero-dimensional systems that 
undergo a noise- induced transition [Q. At small noise 
intensity, the stationary probability density has a single 
maximum at <fi — 0. As the noise intensity crosses the 
critical value a 2 = a/c the stationary probability den- 
sity shows two symmetric maxima about a minimum at 
= 0. This transition involves no symmetry breaking 
since in all cases the average (4> r ) — 0. When sites are 
coupled diffusively, i.e., C = DV 2 , the model ([!]) ex- 
hibits a noise-induced phase transition ||l4|| . For a given 
value of the coupling, as the intensity of the fluctuations 
increases, the systems undergoes a second-order phase 
transition that involves spontaneous symmetry breaking, 
that is, a transition from a state with (<p r ) — to one 
with (4> r ) =/= 0. It is worth noting that when the coupling 
strength parameter D goes to infinity, the location of that 
transition point is the same as in a zero-dimensional en- 
semble, a 2 = a/c fj^j , and that the symmetry-breaking 
in the coupled system can in that case be understood 
in terms of the dynamics associated with the effective 
potential in the zero-dimensional case. 

We show herein that more complex structural changes 
can be obtained by modifying the coupling term C. In 
particular, pattern formation phenomena occur when the 
system develops a morphological instability, that is, when 
a Fourier mode of wavevector k other than k = becomes 
unstable Jl5[ . Drawing parallels with previous literature 
on noise-induced phenomena [^|, [l^, |Tl|, |l2[ , we consider 
a Swift-Hohenberg coupling, that is fllql, 



C = -D (k 2 + V 2 ) 



(4) 



The effect of this coupling can be deduced by applying L 
to a plane wave e lkr , 



Ce ikr = uj(k)e ikr , 



(5) 



where ui(k) = — Z?(/cq — k 2 ) 2 is the (continuous) disper- 
sion relation (we use bold for vectorial quantities and 
italic for their magnitudes). Thus, the largest eigenval- 
ues, Lu(k) — 0, are those associated with Fourier modes 
of wavenumber magnitude k = k$. We will see below 
that this behavior leads to a morphological instability in 
the presence of noise. We stress that a morphological 
instability is a key ingredient of the pattern formation 
mechanism, but the specific functional form of the cou- 
pling term is not. 

In the next two sections we show in detail that our 
model leads to a noise-induced phase transition to pat- 
terned states. We end this section by showing explicitly 
that the transition is not caused by short-time instabil- 
ities H 10, O, The relevant short-time evolution 



equation for the model can be obtained by averaging 
Eq. (|l]) and expanding for small (j> r around the average 
value (6 r ) = 0: 



b r = (—a — a 2 c + C)(j) r . 



(G) 



Writing <f> r (t) as a sum (in an infinite system as an inte- 
gral) over Fourier modes, 



^(t)=J2Mt)e lk - r , 



(7) 



and integrating, we obtain 
(j) r (t) = ^2,4>k exp { [-a - a 2 c + uj(k)] t + ik ■ r} . (8) 

k 

Since u>(k) < 0, all modes decay and there is no short- 
time instability. 

III. MODULATED MEAN FIELD THEORY 

To find the phase diagram for our noisy spatially ex- 
tended system, we introduce a modulated mean field the- 
ory. In order to implement this theory and, in particu- 
lar, to calculate the effect of the coupling operator L in 
the mean field context, we need to explicitly distinguish 
different locations r and r' in a way that requires dis- 
cretization of the system. Since numerical simulations 
also involve discretization, this procedure does not in- 
terfere with the comparisons of theoretical and numeri- 
cal results. With the understanding of the action of the 
translation operator 



(9) 



it is straightforward to deduce the discrete version of the 
Swift-Hohenberg coupling operator, 



C 



-D 



2 

Ax 



sinh 2 

i=l 



Ax_d_ 

2 dxi 



(10) 



where d stands for the spatial dimension, Ax for the 
lattice spacing, and t^- indicates a partial derivative 
with respect the i-component of the position vector r = 
(xi, x 2 , ■ ■ ■ , Xi, . . . , Xd)- As in the continuous case, the 
discrete dispersion relation can be obtained by applying 
the operator diG) to a plane wave e lk r , to obtain 



Lo{k) = -D 



i— 1 



Ax 



(ii) 



Here ki denotes component i of the wave vector k = 
{k\, &2> • • • j ki, . . . , kd). 

Note that as in the continuous problem, uj(k) is non- 
positive for any value of fc, but that in the discrete case 
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it depends not only on the magnitude but also on the 
direction of k. Of particular importance in our subse- 
quent analysis are those modes for which u>(k) = 0. In 
the continuum these are the modes with k = ko, which 
are all those that lie on a continuous hypersurface in re- 
ciprocal space of radius ko around the origin. In the 
discretized system the magnitudes k* of the most unsta- 
ble modes are shifted from fco arL d depend on direction, 
as can be seen by solving Eq. (^Tj). The longest vectors 
such that cu(k*) = lie along the Cartesian directions, 
e.g. (k* ,0,0, ... , 0) and have magnitude 



2 . / k Q Ax 

max k = — — arcsm 

Ax V 2 



(12) 



The shortest lie along a reciprocal space diagonal, e.g. 
■^=(k* , k* , k* , . . . , k*), and have magnitude 



min k* 



2\fd . fk Ax\ 

— — arcsm =■ 

Ax V %Vd J 



(13) 



If koAx < 1 (which it always will be in our analysis) then 
the difference between these two values is smaller than 
3%. It is therefore only a mild approximation to neglect 
the directional dependence of the solutions of u>(k*) = 
and focus on the magnitude, w(fc*) = 0. 

To establish the existence of patterns of a characteris- 
tic length scale, we seek a spatially periodic structure de- 
fined by wavevectors k whose magnitude k is associated 
with the inverse of this length scale. As we shall see, the 
appropriate wavevectors to focus on are precisely those 
of magnitude k* , that is, those for which us(k) = 0. As 
is customary in mean field theories, we make an ansatz 
about the behavior of the field at sites r' other than the 
focus point r which are coupled to it by the operator C, 
one that incorporates an appropriate spatial modulation: 



A{k*) J2 cos \ k ■ (r-r')} 
{k*} 



(14) 



where the sum (or, in an infinite system, the integral) 
is over wavevectors of magnitude k* . Our ansatz thus 
also incorporates the assumption that all modes of this 
magnitude contribute with equal (direction-independent) 
weight A(k*). In the appendix we detail some of the 
steps that show the action of the coupling operator on 
this ansatz, to arrive at the result 



where 



n(fc 



£<j> r = D 1 [n{k*)A(k*) - <j> T ] 



r(d/2 + 1) V 2tt 
is the number of modes of magnitude k* , and 



D 1 = D 



( 2d 



\(Ax) 



k 2 

2 K 



2d 



{Ax) 4 



(15) 



(16) 



(17) 



Substitution in Eq. (|I|) then leads to an equation that 
depends only on a generic site index r that can simply 
be dropped: 



</> = r(0) {-a</> + Di(fc*) [n(k*)A(k*) - (f>}} 

+ [rW0] 1/2 £(t), 



(18) 



and £(t) is zero-centered Gaussian noise (5-correlated in 
time, 



2a 2 



(19) 



Here we have incorporated the fact that the continuum 
delta function S(r — r') has been replaced in the usual 
way by a ratio that contains the Kronecker delta and the 
lattice spacing, S r , r //(Ax) 2 . Henceforth we set Ax = 1. 

The amplitude A(k*) is the mean field quantity that 
must be chosen self-consistently to complete the solution 
of the problem. The stationary probability density for 
the stochastic process (ffa) is 



P (0; A(k*)) = M [n(k*)A{k*)\ (1 + «/> 2 ) 1/2 

Di)(/) 2 — Dxxi(k*)A(k*)<f) 



\ exp <J -^j 



(20) 

where the normalization constant J\f [n(k*)A(k*)] de- 
pends on the amplitude and must therefore be carefully 
included in the self-consistent solution. Self-consistency 
is then embodied in the assumption that n(k*)A(k*) is 
the average value of the field at any point in space, i.e. 
in the requirement that 



n(k*)A(k* 



^(<^;^(fc*))#, 



(21) 



which is appropriate either if A(k*) = and the distri- 
bution is symmetric in <fi, or if n(k*)A(k*) is much larger 
than the (appropriately phased) combined amplitudes of 
all the other modes. The latter occurs if there is an in- 
stability that leads to the formation of a pattern with 
wavenumbers of magnitude k* . 

Since p(0;O) = p(-<f>;0), it follows that A(k*) = is 
always a solution. To find other solutions we expand the 
integral on the right side of Eq. (|2l] ) around A(k*) = 0, 



<j> P {(j);A{k*))d4> = bA{k*) + 0(A i {k*)) , (22) 



where 



dp(<t>;A(k*)) 



dA{k*) 



(23) 



^(fe*)=0 



It follows that self-consistent solutions different from 
A(k*) = are possible, and that the loci that indicate 
the appearance of these solutions satisfy b = n(k*), that 



is, 



^1 P {4>- 0)# = l. 



(24) 
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FIG. 1: Phase diagram of the model obtained from the mod- 
ulated mean field theory. The wide solid line indicates the 
transition loci ( p5[ ) in (D,a 2 ) space for ko = 1 (k* ~ 0.1035), 
a = 1, and c = 3. Inside the ordered region where patterns 
develop we have indicated some contour lines labeled by the 
value of the order parameter S(k*) (see text). Note that as 
the coupling goes to infinity, the transition line tends to the 
value of the noise intensity where the zero-dimensional system 
undergoes a noise-induced transition, a 2 = 1/3. 



The latter condition, which determines the phase dia- 
gram of the model, can also be obtained by geometrical 
arguments H 0. Equation @ can be expressed in the 
following algebraic form: 



2cr 2 cVHJ(-i,0, 



a+Di 



(25) 



where U(x, y, z) is the confluent hypergeometric function 
and K n (x) is the Bessel function of the second kind. 

The derivation is not complete until we confirm, at 
least within the approximations implemented in our 
mean field approach, that the most unstable modes are 
indeed those of magnitude k* , i.e., those for which u(k) = 
0. This is easily ascertained as follows. If we consider an 
ansatz state of the form ( |l4| ) but with a different k, then 
we would arrive at a mean field equation containing uj(k) 
explicitly: 



=T(<p){-(a-uj(k))<p- 



D eff (k) [n(k)A(k) - 0}} 



_ (26) 
where D e //(fc) = D\ + uj(k) (cf. Eq. (A7); here the as- 
sumption has again been made that uj(k) is well approxi- 
mated by w(k). Note that if us(k) = we recover Eq. (18). 
If we now follow the same steps leading to the condition 
( pB| ) which marks the boundary of pattern formation, we 
find that the only modification is that the very first D\ in 
the numerator is now replaced with D e ff(k). The other 
occurrences of D\ are not affected because D\ and a are 
each shifted in such a way that their sum remains un- 
changed. Since D e ff(k) < D\ for all k ^ k* , it follows 



FIG. 2: Order parameter as a function of the noise intensity 
obtained from numerical simulations (squares) for a = 1, c = 
3, D = 5, and ko = 1- The error bars are indicated. The 
dashed line indicates the linear interpolation of the numerical 
data used to locate the transition point, a 2 = 0.5 ± 0.1. The 
points labeled A {a 2 = 0.1), B (a 2 = 2), and C (a 2 = 5) 
correspond to the spatial structures shown in Fig. H. Inset: 
mean field result, with transition point at a 2 — 0.34955 ± 
10" 5 . 



that in order to satisfy Eq. ( |25|) the noise intensity has to 
be greater for other k. In other words, while modes other 
than those of magnitude k* may become unstable, they 
first do so at higher values of the noise intensity. Con- 
versely, for a given noise intensity the coupling must be 
stronger to produce an instability of other wavevectors. 
Thus, as the noise intensity is increased from zero for 
a given coupling, or the coupling is increased from zero 
for a given noise intensity beyond a^, the first and hence 
strongest instability occurs at k* . 

Figure Fy shows the phase diagram of the model in the 
space (a,D) obtained from the numerical solution of 
Eq. ( |25| ) in two dimensions with Ax — 1, a — 1, c = 3, 
and ko = 1, which leads to k* ss 0.1035. Note that as D 
goes to infinity the transition line moves to the value of 
the noise intensity where the zero-dimensional stationary 
probability density becomes bimodal, o 2 c = 1/3. 

Since Eq.(|l|) satisfies the inversion symmetry cf> < — > 
—0, roll-shaped patterns are likely In order to char- 
acterize the structure that develops we use as an order 
parameter the total power spectrum at the most unstable 
modes, 



S(k*) 



(27) 



where, as before, 4>k stands for the Fourier transform of 
the field 4> r and the sum is over all modes of magnitude 
k* . This parameter characterizes the transition from a 
laminar regime (homogeneous) to a convective regime 
(roll-like patterns) that we expect to obtain here. Us- 
ing our self-consistent solution then gives in the modified 
mean field approximation the appropriately normalized 
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FIG. 3: Density plots of the stationary field associated with 
the points A (first panel), B (second), and C (third) of Fig. ^| 
Note that as predicted, for a constant value of the coupling, 
a pattern develops as the intensity of the noise increases. 

relation 

S(k*) = n{k*)A 2 (k*). (28) 

In Fig. [j] we also present several contour lines indicat- 
ing the value of S(k*). Note that for a given value of 
the coupling D, as noise intensity increases, the order 
parameter S(k*) also increases, that is, the stronger the 
noise, the larger is the amplitude of structures associated 
with wavevector magnitude k* . Note also that for suffi- 
ciently large noise intensity, for a given value of a 2 there 
is a non-monotonic behavior of S(k*) as a function of the 
coupling strength, indicating that there exists a value of 
the coupling for which the structures associated with k* 
exhibit maximum coherence. At the same time, the pos- 
sible instability of other modes near k* may affect the 
actual physical appearance of the system, so that these 
effects may not be visually unequivocal. 

Finally, we note that it is common in pattern formation 
discussions to use the so-called flux of convective heat 

J = ^E^ = E^-fc ( 29 ) 

r k 

as an order parameter. If our ansatz state were exact, 
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FIG. 4: Order parameter as a function of the coupling 
strength obtained from numerical simulations (squares) for 
a — 1, c = 3, a 2 = 6, and ko = 1. The error bars are in- 
dicated; the dashed line is a guide for the eye. The points 
labeled A (D — 0, uncoupled system), B (D = f4), and C 
(D = 60) correspond to the spatial structures shown in Fig. ^. 
Inset: mean field result. 

then in the patterned state the two are identical, J = 
S(k*). Differences point to the presence of other unstable 
modes. 



IV. NUMERICAL SIMULATIONS 

In order to check the analytical predictions of the mod- 
ified field theory, we perform numerical simulations of 
Eq. (|l|) on a two dimensional square lattice of 64 x 64 
cells and apply Dirichlet-Neumann boundary conditions 
commonly used in studies of fluids, namely, the field and 
the normal component of the gradient are zero at the 
boundaries. The relevant parameters have the follow- 
ing values in our simulations: mesh size Aa; = 1 that 
gives the system length scale L — NAx = 64, time step 
At = 0.001, ko = 1, which leads to k* as 1.035 and an as- 
pect ratio A = k* L/2tt ~ 10. The other parameters of the 
model are again taken as a — 1 and c = 3. Using these 
values, the onset of bistability of the stationary one-site 
potential occurs at a 2 = a/c = 1/3. We expect the tran- 
sition to pattern formation to occur for noise intensities 
near a 2 if the coupling is sufficiently strong. According to 
the mean-field phase diagram Fig. D = 5 is sufficiently 
large. We use this value of the coupling in subsequent cal- 
culations and figures. Although the initial conditions do 
not matter for the final outcome, the outcome is reached 
more quickly if we start from a target-like pattern of rings 
of width k* . We have ascertained that other initial con- 
ditions, for example a random initial condition, reach the 
stationary state albeit in a longer time. 

We have numerically computed the order parameter 
for several noise intensities, and the results are shown in 
Fig. [2| The order parameter is very small for low noise in- 
tensities until the critical noise intensity is reached, after 
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FIG. 5: Density plots of the stationary field associated with 
the points A (first panel), B (second), and C (third) of Fig. ^ 
Note that, as predicted, for a constant value of the noise, a 
pattern develops as the coupling increases but then becomes 
less distinct as the coupling increases further. 

which the order parameter grows linearly with a 1 . Our 
simulations lead to a critical value (Jq = 0.5 ±0.1. The 
mean field result, also shown in the figure, predicts the 
ordering transition to occur at a 2 = 0.34955 ± 10~ 5 and 
also captures the linear transcritical behavior of S(k*) 
with a 2 , but predicts a slower growth of the order param- 
eter with noise intensity than the numerical simulations. 
Nevertheless, the mean field theory clearly captures the 
full qualitative behavior of the system, especially near 
the transition point. 

We have also computed the flux according to Eq. ( |29| ) 
for all the cases for which we have presented the order 
parameter. We find that both predict essentially the 
same critical noise value for the appearance of patterns, 
and that for our simulation parameters the flux is fairly 
consistently twice as large as the order parameter (be- 
cause of the contributions of other modes to the former). 
This consistency would indicate that even well beyond 
the transition point the most unstable modes dominate 
the flux. 

The spatial pattern that emerges as a 2 increases is 
illustrated in Fig. || using density plots of the stationary 
value of the field for the noise intensities labeled A, B, and 



C in Fig. |2[ a\ — 0.1 (weak noise, no pattern), a 2 B = 2 
(rolls are visible), and cr^ = 5 (strong noise, distinctive 
pattern). 

In order to ascertain the non-monotonic behavior of 
the order parameter with coupling strength at a given 
noise intensity, in Fig. || we present the results for a 
fixed noise intensity a 1 = 6. Simulations are indicated by 
squares, and the inset shows the results predicted by the 
theory. While there are again quantitative differences, 
the qualitative agreement is evident: both curves show 
a clear non-monotonic behavior. The density plots of 
the field associated with points A, B, and C in Fig. || 
confirm this behavior. Of course even as the amplitudes 
of wavevectors of magnitude k* are decreasing, the am- 
plitudes of other nearby wavevectors may be growing. 
Therefore, while the lightening of pattern C on the gray 
scale reflects the decrease in the order parameter, a vi- 
sual perception of some loss of distinctness could be due 
to an admixture of other wavevectors. 



V. CONCLUSIONS 

We have shown by means of a modulated mean field 
approximation and numerical simulations that the mech- 
anism for noise-induced phase transitions introduced by 
Ibanes et al. |l4|] can be extended to pattern formation 
phenomena. In contrast with previous work on noise- 
induced patterns, no short-time instability is required 
to generate these spatial structures. As a consequence, 
this transition is independent of the noise interpreta- 
tion, as has been shown for noise-induced phase tran- 
sitions |L8|]. For example, for the ltd interpretation we 
have checked numerically that this model exhibits the 
same noise-induced patterns, but at a lower critical noise 
intensity. Furthermore, no rc-cntrancc phenomenon oc- 
curs as the fluctuations grow in intensity. Indeed, in our 
system stronger noise leads to increasingly ordered struc- 
tures. The mean field theory has allowed us to charac- 
terize the model by computing its phase diagram. We 
have also performed numerical simulations that confirm 
the qualitative validity of the theoretical analysis. We 
find that as the intensity of the fluctuations increases, a 
rotationally symmetric roll-shaped pattern appears. The 
pattern is characterized by the most unstable modes of 
the system, those with wavevector of a magnitude k = k* 
that is explicitly predicted by the mean field analysis. 
Both the theoretical and numerical analyses show that as 
the coupling between sites goes to infinity, the transition 
to pattern formation occurs at the same point where the 
zero-dimensional system presents a noise-induced tran- 
sition. We have also shown that for sufficiently strong 
noise intensity the order parameter for the system is non- 
monotonic as a function of the coupling strength. Thus, 
for sufficiently strong noise there exists an optimal value 
of the coupling such that the patterns of characteristic 
size 2n/k* are maximally coherent. 
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APPENDIX A: IMPLEMENTATION OF THE 
MODIFIED MEAN FIELD THEORY 



We begin by exhibiting in some detail the dependences 
associated with the ansatz field ([l4]) and the action of C 
on it. For example, for an r' that is m lattice sites away 
from r = (x\,X2, ■ ■ ■ ,Xd) in the direction j the ansatz 
reads 



(x%,X2, . . . , xj+mAx, 



,Xd) = 53 A(k*) cos(mAxkj). 
{k>} 

(Al) 

For an r' that is in the immediate positive diagonal lo- 
cation away from r we have 

4>{x\ + Ax,X2 + Ax, . . . , Xj + Ax, . . . ,Xd + Ax) 

= 53 A ( k *) cos[Aa;(fci + k 2 H h k d )}. 

{k*} 

_ (A2) 
Next, to apply the discrete version <M) of £ we must elu- 



cidate the effect of the operators 

on the field 4> r for n = 1,2. With n = 1, we use the re- 
lation 2sinh 2 (j//2) = [cosh(y) — 1] and note that 



53 cosh ( Ax 



<t>r 



1 



dxi J 2 

Xi — Ax, X2, ■ ■ ■ ,Xj, 



[(f)(xi + Ax, X2, ■ ■ ■ , Xj 



{Xi,X 2 , ■ 

{xi,x 2 , ■ 
(xi,x 2 , ■ 
l )(xi,x 2 , 



Xd) 

Ax, ... ,x d ) 
Ax, ... ,x d ) 
.,x d + Ax) 
..,x d - Ax)} 



(A3) 



By using Eq. (Al) in this last equation we obtain 



d 

£>sh(A^ 

i— 1 



53 A(k*)J2 C0S ( k * A x)- 

{k*} i=l 



2 (A4) 

As for n — 2, we note that 4sinh (y/2)sinh (z/2) — 
[cosh(y) — 1] [cosh(z) — 1] and, in turn, cosh(y) cosh(z) = 
i[cosh(j/ + z)+cosh(y — z)]. The latter combination leads 
to contributions that involve both forward and backward 
translations in different spatial directions. This is easily 



visualized by noting explicitly that 

2 

d 



^ cosh (Ax— 



5>sh(Ax( ^- + ^- 



cosh ( Ax 



_d_ 

dxi 



dxi dx^ 
d 

dXj 



(A5) 



Notice that for the d cases where with i = j, the second 
term on the right hand side leaves the field at the original 
site r. The field at the original site is not represented by 
the ansatz assumption, and therefore we must subtract 
the d "spurious" terms produced by the ansatz state and 
add d times the field cf> r . This procedure leads to, 



53 cosn ( A 



d 

dxi 



53 A(k*) 



53 cos(kiAx) 



\i=l 



(A6) 

Note that we have taken advantage of the directional 
insensitivity of k* . 

Use of Eqs. ([X|) and (|A6|) in Eq. @ then leads to 
the following approximation for the term containing the 
SH coupling operator: 



C4> r =D X 53 Ak*)- 

\{k*} 

=Dr I 53 A(k*) - 

\{k*} 



53 A{k*Mk*) 
{k»} 



(A7) 



where the last term vanishes because tu(k* 

v 2 



0, and 

where £>i = D [(^ - 

Finally, since the summand in Eq. (A7) is indepen- 
dent of the direction of the k* , the sums simply give the 
number of terms in the sum (or the appropriate integral 
form) times the summand. Simulations always involve 
a finite system of N d sites (i.e., of volume (NAx) d ), so 
that the allowed modes themselves form a discrete set, 
with each component separated from the next one by 
an interval Ak — tJj-- One way to count the number 
of modes n(k*) in the sum is to construct a ring of ra- 
dius min k* (which we shall simply call k* following the 
discussion surrounding Eqs. (|l^) and (|l3|)) of thickness 
Afc = 2n/N Ax, and to consider all the modes that lie in 
this ring. We can then estimate this number by calcu- 
lating the number of cells of volume (2tt / N Ax) d in the 
ring: 



n(fc*) 



dir d / 2 (Nk* 



r(d/2 + l) V 2tt 



d-i 



(A8) 



8 



Although variations in the particular way of counting are 
possible, for sufficiently large N the differences are small. 
Thus we finally arrive at the mean field approximation 

£<p r = D 1 [n(k*)A(k*)-(j> r }. (A9) 
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